Front Propagation in Chaotic and Noisy Reaction-Diffusion Systems: 

a Discrete-Time Map Approach 
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We study the front propagation in Reaction-Diffusion systems whose reaction dynamics exhibits 
an unstable fixed point and chaotic or noisy behaviour. We have examined the influence of chaos and 
noise on the front propagation speed and on the wandering of the front around its average position. 
Assuming that the reaction term acts periodicaUy in an impulsive way, the dynamical evolution 
of the system can be written as the convolution between a spatial propagator and a discrete-time 
map acting locally. This approach allows us to perform accurate numerical analysis. They reveal 
that in the pulled regime the front speed is basically determined by the shape of the map around 
the unstable fixed point, while its chaotic or noisy features play a marginal role. In contrast, in 
the pushed regime the presence of chaos or noise is more relevant. In particular the front speed 
decreases when the degree of chaoticity is increased, but it is not straightforward to derive a direct 
connection between the chaotic properties {e.g. the Lyapunov exponent) and the behaviour of the 
front. As for the fluctuations of the front position, we observe for the noisy maps that the associated 
mean square displacement grows in time as f^' ^ in the pushed case and as t^'* in the pulled one, in 
agreement with recent findings obtained for continuous models with multiplicative noise. Moreover 
we show that the same quantity saturates when a chaotic deterministic dynamics is considered for 
both pushed and pulled regimes. 

Pacs numbers: 05.45.-a, 05.45. Ra, 47.20. Ky, 68.10.Gw 



I. INTRODUCTION 

In the last years the study of front dynamics has gained 
more and more relevance in many different fields, such as 
chemistry [Q, biology [||, combustion |^. In physics Q 
the problem of front propagation is generally related to 
Reaction-Diffusion (RD) systems, and to the identifica- 
tion of the different dynamical regimes present in the 
model under study. 

In this case the focus is often on the situation where RD 
fronts connect a stable state to an unstable one. Consider 
for instance the prototype equation for front propagation 



wiin = mini;(Ai) = 2/Df(0) 



and 



2i L> max 

V 0<«<1 u 



(2) 



(3) 



dtu — Ddr 



f{u) 



(1) 



where u — u{x, t) , D is the diffusion coefficient and / is a 
continuous function with two fixed points, u = 1, stable, 
and u = 0, unstable. This equation is usually referred 
to as Fisher-Kolmogorov-Petrovsky-Piskunov (FKPP) 
Equation [^,||. 

It is well known that if one considers an initial con- 
dition m(x, 0), different from zero in a bounded spatial 
region, a front develops connecting the two fixed points. 
The asymptotic leading edge, i.e. the semi-infinite region 
ahead of the front, has typically an exponential shape of 
the type exp(— /ix + Xt). In principle, depending on the 
value of the spatial decay-rate ^, the front speed can take 
a continuous set of values, namely v{iJ,) = X/fi 0. 

However, for sufficiently localized initial conditions a 
unique speed V{ is selected, which is always bounded in 
the range [wii„ , Wmax] , with 



If the function f{u) is concave, it is possible to prove 
that the propagation speed of the front coincides with 
the linear one, V{ = vun P|- In literature this case is 
often referred to as pulled dynamics. The nonlinearities 
present in the system are not dynamically relevant, and a 
linearization about the unstable state is enough to fully 
describe the propagation of the front. As a result the 
front is indeed "pulled" by the spreading and growth of 
linear perturbations in the leading edge ||8| . On the other 
hand, the function f{u) being not concave is a necessary 
condition to have V{ > wnn f^. This in turn corresponds 
to the pushed regime, where the nonlinearities are rele- 
vant to the dynamics and the front is "pushed" by its 
internal part ||^. 

Such a dynamical behaviour can also be studied when 
inhomogeneous environments are considered. A typi- 
cal example is the photosensitive Belousov-Zhabotinsky 
chemical reaction. This system allows for an experimen- 
tal realization of external noise, achieved via fluctuating 
illumination conditions [nOJ. The random environment 
can be mimicked by multiplicative noise terms in the 



corresponding RD equation - see for example [[11|^4[. 
Under the effect of noise, the motion of the front can 
be decomposed into a global drift, characterized by a 
noise-renormalized speed, plus fluctuations, affecting the 
position of the front. Particularly the wandering of the 
front around its average position has been proven to be 
diffusive in the pushed case p2| and subdiffusive in the 
pulled one |14| , with an associated mean square displace- 
ment growing in time respectively like t^/^ and t^/^. 

Another natural way for erratic behaviors to occur in 
RD systems is via a chaotic underlying process. Even 
though so far only a limited number of papers has been 
devoted to the influence of chaotic dynamics on front 
propagation jT^.Q, chaos is extremely relevant in RD 
systems describing chemical reactions 0, as well as in 
more general pattern forming systems [p5[ . Moreover, in 
the context of spatially extended chaotic systems many 
concepts have been borrowed from the study of front 
propagation into unstable steady states to describe the 
spreading of information ijl^l . 

In this paper we aim at taking up this issue again, by 
studying the front propagation in systems where an un- 
stable fixed point is present and in addition the reaction 
dynamics is chaotic (or noisy). In the chaotic case the 
role of the stable fixed point in the usual FKPP problem 
is played by the chaotic phase. The most natural way 
of realizing this situation is to consider the case u € i?'' 
with d > 3 in the FKPP Equation, that is 



dtu = DdxxVi + f (u) 



(4) 



Here u = (ui, . . . , Ud) and f = (/i, . . . , /d) is such that 
u = is an unstable fixed point. Then, omitting diffu- 
sion, equation du/dt = f (u) can be chaotic and therefore 
we expect front solutions connecting the unstable state 
u = to the chaotic state to be realizable. 

However, as it is well known, there is an easiest way 
of introducing chaos in Eq. d^), without resting on the 
generalization to vectorial fields. Since our interest is 
in the qualitative effect of chaos in the reaction terms 
of (|l|), we can reduce ourselves to studying the simplest 
chaotic system, i.e. Id discrete-time maps. As we shall 
see. Id discrete-time maps allow for the analysis of both 
chaotic and non-chaotic dynamics, either in the pushed 
or in the pulled propagation regimes. Also, they can 
easily be generalized to include noise in the system, per- 
mitting thereby a comparison between deterministic and 
noisy dynamics. 

Detailed numerical simulations of such systems show 
that in the pulled regime the front speed is basically de- 
termined by the unstable fixed point, while the chaotic 
(noisy) character of the reaction dynamics plays a lim- 
ited role. In contrast, in the pushed regime the pres- 
ence of chaos (noise) has some relevance, even if the rela- 
tionship between chaotic properties and the front speed 
appears far from being trivial. The differences between 
the chaotic and the noisy situations are much more ev- 
ident in the wandering of the front around its average 



position. The fluctuations of the front induced by the 
chaotic dynamics appear to be bounded on the examined 
time scale, while the presence of noise induces diffusive 
or sub-diffusive behaviour, in the pushed and pulled case 
respectively. 

The discrete-time map approach is introduced in Sec- 
tion II, while the specific map models under study are 
introduced in Section HI. The numerical results concern- 
ing front speed and front fluctuations are discussed in 
Section IV, and flnally some concluding remarks are pre- 
sented in Section V. The details on the integration al- 
gorithm and the estimation of the speed bounds for the 
discrete-time map approach are reported in Appendix A 
and B, respectively. 



II. DISCRETE-TIME MAP APPROACH 

With respect to Eq. (Q), there exists a simple way for 
introducing discrete-time maps. This rests on consider- 
ing the case when the reaction term acts in an impulsive 
way. 



/(m) = Atg{u) Y^ 6{t - nAt) 



(5) 



where g{u) is a function whose precise form is not rele- 
vant at this stage, and At is the period of the forcing. 

Consider now a kick given at time t. By direct integra- 
tion of Eq. (|l]) between t and t + e, in the limit e ^ 0, 
we obtain 



where 



u{x,t + 0+) = F(u(a;,t)) 



F{w) =w + Atg{w) 



(6) 



(7) 



F{v) is the reacting map. In other words we have written 
the evolution of the field between two successive kicks in 
terms of a simple map. 

The remaining evolution of the field can be calculated 
noticing that from one kick to the next one, the evolu- 
tion of the field is "free" because the impulsive reactive 
term by definition does not act. Indeed between t -|- 0+ 
and t + At, equation (m reduces to the diffusion equa- 
tion dtu ~ DdxxU, whose solution is known and allows 
for writing the complete evolution as: 



u{x,t + At) 



1 



V4:TTDAt 

1 

V^TTDAt 



dy e 'toAtu(x — y,t + 0^) 

dy e '^o^F[u{x ~ y,t)] . (8) 



Let us notice that the above equation is nothing but the 
discrete-time version of the well known Feynman-Kac for- 
mula |19|| . A similar approach has been recently intro- 
duced in pO| | to study the two-dimensional dynamics of a 
front separating a chaotic state from a stable steady one. 



We are only left with the evaluation of the convolution 
integral (^. Its numerical estimation can be performed 
employing a quite efficient and accurate algorithm re- 
cently introduced in [gl|. A detailed description of the 
algorithm is given in Appendix A. Once the map is given, 
then its successive iterations account for the whole time 
evolution of the field according to Eq. (||). We remark 
that Eq. (j^ is exact, i.e. it is not an approximation for 
small At if Eq. (ph holds. For the sake of simplicity in 
the following we shall set At = 1. 

Of course the evolution of the field will be chaotic or 
not depending on the map. It is well known that the 
appropriate dynamical indicator to discriminate between 
chaotic and non-chaotic dynamics is the maximal Lya- 
punov exponent A, which characterizes the divergence in 
time of nearby orbits. This exponent is estimated by ap- 
plying a standard evaluation scheme [g2| to the evolution 
of an infinitesimal perturbation 6u{x,t) of the reference 
orbit in the tangent space 



Su{x, t + 1) 



1 



xF'[u{x 



dy e 4D 

y,t)]Su{x~y,t) 



(9) 



The maximal Lyapunov exponent A is then defined as 

(10) 



A=limilnf^ 
t^oot \\Su{x,0)\\ 



where ||(Su(a;, t)|p = / \du{x,t)\'^dx, A > representing a 
signature of chaos. 



where c — ui{a ^ (3) , d = 1 — j and U2 — {d — c) / (P — 7) 
(see Fig. ||a). 

The fixed point uo = 1 is stable provided that I7I is 
smaller than one. As we shall s ee, for this system the 
linear speed is uun = 2\/D\na while the slop e of the 
exponential part of the leading edge is /xun = ^Ina/L). 
It is clear that we are always in the pulled situation if 
a > P, while if f3 is bigger than a pushed situations can 
be observed. 

As a second reference model for non-chaotic dynam- 
ics we also introduce Map B. Now we set U2 = 1/2 and 
F{u2) = 1, while c = 1-/3/2, ui = c/(a-/3), d = I-7/2 
and uq — d/{l — 7). Also in this case we choose I7I < 1 
in order to have a stable fixed point in uq. This map is 
shown in Fig. ||&. The condition a < j3 allows again for 
pushed situations to be in principle observed. 




u, u u, Uj u„ u 

FIG. 1. Panel a): Map A. Since a > [3, the expected front 
propagation corresponds to the pulled case - see text. Panel 
b): Map B. In this case a < /3 makes it possible to observe 
pushed dynamics. 



III. THE MODELS 

We describe first the two deterministic non-chaotic 
maps, which turn out to be useful to illustrate the general 
approach and to identify the different regimes present 
in the system. Starting from those models two general- 
ization will be proposed, in order to assess whether the 
effects of noise or chaos on the system are relevant. 



B. Noisy Models 

In order to mimic a random environment, in |12|-|1^ 
the following noisy RD system was considered: 



dtu = DdxxU + f{u) + g{u)ri 



(12) 



with 



A. Deterministic Non Chaotic Models 



f{u) = u{\ — u){a + m), g{u) — u(l — u) 



(13) 



The first map we introduce has the purpose of repro- 
ducing the usual FKPP behaviour. As mentioned above, 
this is given by equation (Q) where the function / is cho- 
sen in such a way that u = 1 is a stable fixed point and 
M = is an unstable one. Inspired by this feature, we 
propose a deterministic non-chaotic map with an unsta- 
ble fixed point in u = [i.e. F(0) = and |F'(0)| > 1), 
and a stable one uq, {i.e. F{uq) = uq with |F'(wo)| < 1). 
More specifically we define Map A as 



F{u) 




< u < ui 

Ui < U < U2 
U2 < U < 1 



(11) 



The parameter a is a control parameter, by tuning which 
it is possible to change the stability properties of the in- 
vaded state, thereby exploring both pushed and pulled 
dynamics [121. The noise ry = r]{x,t) is Gaussian, spa- 
tially and temporally (5-correlated, and because of the 
chosen g{u), it does not affect the fixed points and 
1. The multiplicative noise term present in ( |l2| ) was 
first introduced phenomenologically as external noise in 
|12| , [l4| , and then rederived in a broader context in p4[ . 
Specifically, in Il2,n4|, due to the continuous nature of 
the model (E2) it was proven that the appropriate pre- 
scription for the evaluation of the noise term was the 
Stratonovich one |23|. 



Relevant questions about such systems are typically 
related on one side to the computation of the renormal- 
ization due to noise of the propagation speed of the front, 
and on the other side to the identification of the general 
properties of the wandering of the front around its aver- 
age position. 

This scenario rests on the decomposition of the prop- 
agation of the front, solution of (112), into a global drift 
plus fluctuations. The global drift is associated to an av- 
erage deterministic front, where the average is taken over 
the different realizations of the noise. This front obeys 
a deterministic field equation which can be obtained 
through a standard procedure from equation ( n2 l ) (s ee for 
example [^). In the simple case of the choice (|llt), this 
deterministic equation reduces to an equation with the 
same reaction term of (03) but with noise-renormalized 
parameters, implying thereby a renormalization of the 
front propagation speed ||ll|. Notice that this effect is 
related to the assumed Stratonovich prescription in the 
evaluation of the multiplicative noise term. 

Of course the single realizations of the noise induce 
fluctuations in the shape of the front, which manifest 
themselves as fluctuations of the speed around its aver- 
age value. This process induces in turn a wandering of 
the front around its average position. This can be charac- 
terized by the root mean square displacement of the front 
position A(t). As a result, the pulled and pushed regime 
differ noticeably. In particular, in the pushed regime the 
usual diffusive behavior is proven [OUO] i.e. A(t) ^ t^l"^ 
while in the pulled case a theoretical and numerical anal- 
ysis shows that subdiffusion occurs, with A(i) ^ t^/"' 

0- . . . 

Having in mind all this, and starting from the deter- 
ministic non-chaotic models introduced in the previous 
subsection, namely Map A and Map B, we introduce 
the corresponding noisy maps. Even if a direct compari- 
son between the continuous model and our discrete-time 
maps is beyond the scope of the present paper, never- 
theless we aim at assessing if the main features of the 
continuous FKPP equation with multiplicative noise are 
still observable for periodically forced discrete systems. 

As a first choice, we consider the stochastic version of 
Map A reported in Eq. (O). We shall refer to this model 
as Map NA and we expect that it will be qualitatively 
equivalent to the noisy FKPP Equation, Eq. (|l3|). Of 
course the insertion of noise can be performed on any of 
the three segments of the map, identified by their slopes 
a, /3, 7. Our strategy will be of inserting noise on one 
of them, by keeping fixed the other ones. For instance 
we can consider a situation where both (3 and 7 are fixed 
while a is chosen randomly at any space-time point as 



a{x, t) = uq + r]{x, t) 



(14) 



Here the noise 77 is distributed according to a fiat distri- 
bution defined in the interval [— v4, A], and is (5-corrclated 
both in space and time. The model is constructed in such 
a way to maintain u — Q and u = 1 as fixed points. 



The second type of noisy map that we analyze can be 
introduced starting from Map B. In this case the noise is 
added to the parameter 7 and it has again the same flat 
bounded distribution of Map NA and the same correla- 
tor. In this case, the restriction that I70 ± A\ be smaller 
then one forces the map to have a unique unstable flxed 
point u = but an infinite number of different stable 
fixed points u = uo(7)- As a matter of fact the inte- 
rior region of the front, developing from the state u — Q, 
oscillates stochastically around the average flxed point 
(uo) — d/{\ — 7o). We shall refer to this model as Map 
NB. It will be a benchmark to understand what the ef- 
fects are of the noise, when it affects the saturated part 
of the front. 

C. Chaotic Models 



Modifying Map B it is possible to introduce a chaotic 
version of the FKPP Equation, which we shall term Map 
CB. This map is deflned by keeping it identical to Map B 
for M < M2 = 1/2 and modifying the part in the interval 
1/2 < u < 1 in the following way: 



Flu] 



ju + d, 
Su + e, 



U2 < u < U3 
U3 < u < 1 



(15) 



where 7 = 1/(1/2 ~ U3) < 0, S ^ 1/(1 - M3) > 0, 
d — I — 7/2 and e = 1 — 6, while U3 is arbitrar- 
ily chosen with the restriction to belong to the interval 
1/2 < U3 < 1. Map CB is represented in Fig. 0. 



F(u) 




Ui UjUo U3U 

FIG. 2. Map CB. A variation of the value of us induces a 
modification of the corresponding Lyapunov exponent. 

Notice that the flxed points of the map (w = and 
Mo) are both unstable, because now the value of \a\ and 
I7I are bigger than one. With these conditions, the map 
appears to be always chaotic, i.e. the maximal Lyapunov 
exponent A is positive. 

As shown in the next section, in this case a localized 
initial disturbance of the unstable state u = will grow 
and spread along the system. As a consequence, a front 
will develop connecting the unstable flxed point m = to 
a chaotic phase, playing the role of the invading state. 



IV. NUMERICAL RESULTS AND DISCUSSION 

Let us first present some qualitative results concern- 
ing the main features of deterministic non-chaotic fronts 
{Map B), noisy fronts {Map NB), and chaotic fronts {Map 
CB), in both the pushed and pulled regimes. 

In all the examined cases, the chain length was L = 
30, 000, and it was initially set everywhere to 0, with ex- 
ception of a few sites in the center of the chain, which 
were initialized with a disturbance amplitude of 0(1). 
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FIG. 3. Front evolution for the FKPP equation with Map 
B in a) the pushed regime and h) the pulled one. In the 
pushed regime the values of the parameters are a — 1.001 
and l3 = 5.0, while in the pulled one we set a = 1.8 and 
(3 = 2.5. In both cases D = 4 and 7 = —1/2. 

In Fig. 0a and pp we show the deterministic non- 
chaotic behaviour of pushed and pulled fronts respec- 
tively. The corresponding dynamics is selected by chang- 
ing the relative weight of the parameters a and /3, as 
explained in Section III A . The realization of the pushed 



or pulled regime is checked by direct measurement of the 
front speed (see the precise definition in the next Sec- 
tion). The linear speed corresponding to case shown 
in Fig. gp is uun = 0.126, while the measured value 
is V{ = 1.444, confirming thereby the realization of the 
pushed regime. As for Fig. ^6, we measure V{ = vun, with 
wiin = 3.066, corresponding indeed to pulled dynamics. 
In the pushed case it is evident an abrupt jump from 



an initial evolution where the dynamics is ruled by the 
linear mechanisms (characterized by a Gaussian shaped 
perturbation propagating with velocity vnn) to a situa- 
tion where nonlinearities set in and the front saturates 
in its central part and starts propagating with a velocity 
V{ > wiin- In contrast, in the pulled situation (depicted 
in Fig. 0&) the transition from the initial Gaussian per- 
turbation to a saturated propagating front is smoother, 
since now the only mechanism responsible for propaga- 
tion is the linear one and the speed does not increase 
above the linear value. 

Noisy front propagations as given by Map NB are 
shown in Fig. Ba and 06 for pushed and pulled dynamics 
respectively. In this case the values of the parameters are 
the same as in Fig. 0, with the only exception of the 7 
parameter, which is the one affected by the noise. 
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FIG. 4. Front evolution for the FKPP equation with Map 
NB in a) the pushed regime and h) the pulled one. In the 
pushed regime the values of the parameters are _D = 4, 
a = 1.001 and /3 = 5.0, while in the pulled one we set D = A, 
a = 1.8 and /3 = 2.5. The noise is set on the parameter 7 
with A = 0.45, with 70 = -1/2. 

It is remarkable the fact that there seems to be no 
effect of noise on the value of the speed. As a matter 
of fact, the measured speeds have the values vi = 1.444 
(against the linear value wun = 0.126) in the pushed case 
and again Vf = vun = 3.066 in the pulled one. These 
values coincide with the ones obtained in the determin- 



istic non-chaotic dynamics. However, by analogy with 
the continuous model, we could naively expect a renor- 
malization of the front speed to occur in both pulled and 
pushed regimes. The fact that this does not seem to 
be the case is related to the intrinsic discreteness of our 
models. In the next Section we shall analyze this point 
in greater detail. 

Our results for the chaotic models are shown in Fig. H, 
again for both pushed and pulled dynamics. 
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FIG. 5. Front evolution for the FKPP equation with Map 
CB with U3 = 1 in a) the pushed regime and b) the pulled one. 
The parameters q, /3,and D are again chosen as in the pre- 
vious figures, while the parameter 7 is now set to the value 
7 = — 2. The snapshot shown in Panel 6 is taken at time 
i = 90. 

As it is evident by comparison of Fig. ^ and ||, chaos 
and noise seem to affect in a qualitatively equivalent way 
the front structure, at least when the noise is set on the 
7 parameter. This could suggest the possibility of easily 
building up an effective noisy model for a general under- 
lying chaotic dynamics. However as we shall see in the 
next Section, an effective equivalence between the two 
models is far from being trivial, due to the different role 
played by noise and chaos on the renormalization proper- 
ties of the speed of the front. This is clear in the pushed 
regime, where the front speed takes the value vi = 1.276, 
smaller than in the corresponding non-chaotic and noisy 
models, indicating a deep difference between the respec- 



tive dynamics. 

Finally notice that all the fronts shown in the pre- 
vious figures manifest the same dynamical evolution at 
short times up to t* — C'(l/ln(a)). As a matter of fact, 
in the pushed case the front evolution associated to the 
three maps become distinguishable at times of the or- 
der oi t — 4,400, with a — 1.001, while for the fronts 
in the pulled regime the evolution of the three maps are 
practically identical up to time t = 12, since now the a 
parameter is much bigger (namely, 1.8). Afterwards it 
becomes possible to distinguish the different evolutions 
in both situations. 

In order to investigate at a more quantitative level 
the different dynamical behaviour of fronts in presence 
of noisy or chaotic dynamics, we analyze now the front 
propagation speed and the fluctuations of the front posi- 
tion. 



A. Speed of the Front 

It is not difficult to show that for the discrete-time map 
version of the FKPP Equation, the propagation speed Vf 



is always bounded in the interval [vy, 



111: ^inaxj 7 



where 



V\u 



= 2^DMF'iO)) 



and 



2Wi:'ln 



F{u) 
max 

0<M<1 u 



(16) 



(17) 



Justification of (Hq) and (Ww is given in Appendix B. 
As mentioned before, if Uf = uun, then the dynamics is 
pulled, while if V{ > vn^ the corresponding dynamics will 
be pushed. 

From the numerical point of view, the measurement of 
the speed has been performed in the following way. After 
having initialized to zero all the chain apart a localized 
disturbance u{x, 0) = C(l) for — ^ < 2: < <^, at each time 
step the rightmost r{t) and the leftmost position l{t) of 
the front are considered, 



r(t) — maxja; | u{x, t) > ■&} 
l{t) = min{x I u{x,t) > d} 



(18) 
(19) 



where ?? is a preassigned threshold. The position of the 
front is simply a:j(t) = r{t) — ^ or x^(i) = ^ — l{t), de- 
pending on which of the two moving fronts is considered, 
and the speed results in 



Vi 



lim 



x^f{t)-x\{t) 



2i 



(20) 



We have verified that vi does not depend on the chosen 
threshold and obviously the time limit reported in ( |2(]| ) 
should be taken after the limit L — > 00. 

We present now a more detailed investigation of the 
dependence of the speed on both noise and chaos. Con- 
sider first Fig. p. Our main result is that the "structural" 



features of the map (its concave or convex character) are 
sufficient to provide information about the type of front 
propagation observable and the chaotic, noisy or non- 
chaotic nature of the dynamics seems to be secondary. 
This conclusion relies on the measurement of the front 
speed for Map B, Map NB, and Map CB with U3 = 1. 
We have performed the simulation by keeping (3 — 2.5 
and 7 — —1/2 fixed, and changing the a values in the 
interval 1 < a < 2. As shown in Fig. p^ Uf > v\\a for 
a < 1.4, and therefore for a > 1.4 the linear mechanism 
prevails on the nonlinear one and the front is pulled from 
the linear instability of the leading edge. The same choice 
of parameters has been done also for maps CB {u^ — 1) 
and NB in such a way that the three maps essentially 
coincide for u < M2 = 1/2. 




FIG. 6. Front speeds vi as a function of the parameter a 
for the Map CB (us = 1) (triangles), B (circles) and NB (di- 
amonds) with j3 = 2.5 and D = 4. The dashed line indicates 
Umax = 2VD\n2 and the solid one vun = 2VDlna. The 
length of the chain is L = 10, 000, the front has been followed 
for a time t = 1, 600, and the average is performed over 400 
different initial conditions. 

As we can observe from Fig. || the speeds of the 
fronts are almost identical in the three cases. Only in 
the strongly pushed regime, that is for values of a very 
close to 1, the propagation for the chaotic model appears 
to be slightly slower than in the other models, but the dis- 
crepancy is indeed very small. For values of a as small as 
a = 1.001 the observed values for the velocities obtained 
for both Map B and NB are vi = 1.965, while for Map 
CB (with U3 = 1) we obtain Wf = 1.884, corresponding 
to a discrepancy of the order of 4%. 

In order to better analyze these findings, we investigate 
the dependence of the speed separately as a function of 
the intensity of the noise A and of the Lyapunov coeffi- 
cient A. 

As far as the effect of noise is concerned, we consider 
first Map NB. We set noise on the three parameters of 
the map, i.e. a, (3 and 7, but the corresponding changes 
of the speed appear basically negligible in all cases. 

More interesting is the case corresponding to Map NA. 



In this situation, there is still no effect of noise in the 
pulled regime, while in the pushed case we find a decrease 
of the speed with the noise amplitude, for noise added 
to any of the three parameters. Anyway, the strongest 
decrease has been observed when the noise affects the 
parameter 7. This case is the one reported in Fig. 0. At 
this point it is worth to remark that the very weak veloc- 
ity renormalization due to noise observed in the present 
context is not in contradiction with the results reported 
in |l|], where a noticeable variation of the front speed 
with the noise amplitude was found. 



This is so because the intrinsic discrete nature of our 
system does not leave room to any ambiguity in the defi- 
nition of the multiplicative noise term (the so-called Ito- 
Stratonovich dilemma |23| ) . As a matter of fact, in the 
limit At —* our model would correspond to a contin- 
uous model with multiplicative noise inserted according 
to the Ito's prescription. 




3.05 



FIG. 7. Dependence of the front speed V{ on the intensity 
of noise A in Map NB. Here the noise is inserted on the param- 
eter 7, with 70 — —1/2. The values of the other parameters 
are a = 1.001 and /3 = 5 (strongly pushed), and D — 4. 



The second point worth of investigation is the effect 
of chaos. In this case we performed a measurement of 
the speed of the front as a function of the Lyapunov co- 
efficient in the strongly pushed regime. A different de- 
gree of chaoticity is obtained by tuning the value of U3 
in Map CB, and can be evaluated via the corresponding 
Lyapunov exponent A. Our results are plotted in Fig. @. 



The speed exhibits a strong decrease as a function of 
A, which seems to indicate a behaviour similar to noisy 
case. However, by direct comparison of Figs. ^ and g, it 
appears evident that chaos affects the system in a quite 
more relevant way then noise. 




FIG. 8. Front speeds Vf as a function of A for the Map 
CB. The values of the parameters are in this case D = A, 
a = 1.001, /3 = 5, Ml = 0.3751. The length of the chain is 
fixed for all the simulations to L = 30, 000, the front has been 
followed for a time i = 5, 000, and the average are performed 
typically over 1,000 different initial conditions. 



B. Fluctuations of the front 



Now we consider the root mean square displacement 
of the front position Xf{t), 



A{t) = ^{x}{t)) 



'(t))' 



(21) 



where the average is taken over different initial condi- 
tions for the chaotic case and over many distinct noise 
realizations for the noisy case. Different scalings for A(t) 
are observed depending on the type of map and on how 
the noise enters the dynamics. 

Let us examine first Map NA. We have studied the 
three different cases corresponding on applying noise on 
a, P or 7. In the pulled case the noise does not induce 
any wandering of the front when it is added to /3 or 7: 
In such cases A(t) — > const, in the limit t —^ 00. In 
contrast, if the noise is applied to the a parameter, sub- 
diffusive behavior is indeed observed. These results can 
be justified recalling that in the pulled case the dynamics 
of the front is determined in the leading edge. This re- 
gion corresponds to small u values, and therefore only a 
(stochastic) change in a is expected to affect the behavior 
of the system. 
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FIG. 9. Fluctuations A(t) as a function of time in a log-log 
scale for Map NA, with noise on the a-parameter. Namely 
ao = W, A = A, f3 = 1.0, 7 = -0.2, and D = 4. The data 
have been averaged over 4, 400 different noise realizations for a 
chain of length L = 60, 000, with bounded noise. The dashed 
line refers to a power law t'' with exponent q — 1/4. In this 
case the linear speed is vun = 6.0697, while the measured one 
is Vf ~ 6.069: we are clearly in the pulled case. 

On the other hand, in the pushed regime we observe 
that in all three cases the noise has the same effect on 
the front wandering and it leads to a diffusive behaviour 
for the front positions (see Fig. |I^. This effect is also 
reasonable since now the front propagation is also related 
to the regions where the field takes values of 0(1), and 
not only to the leading edge. Therefore the effect of the 
noise will be equally relevant, independently of what pa- 
rameter it is applied to. 
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FIG. 10. Fluctuations A(f) as a function of t in a log-log 
scale for Map NA, with noise added to the parameter a (with 
values of the parameters uq = 2, A = 1, /3 = 6, and 7 — —0.2), 
l3 (with values of the parameters a — 2, j3o — 6, A = 3 
and 7 = —0.2) or 7 (with values of the parameters a = 1.1, 
(3 = 2.5, 70 = -0.5, and A = 0.4); always D = A. The solid 
line represents the power law f' with q = 1/2. We are in the 
pushed situation since V{ > iiun for all three examined cases. 
Namely for noise acting on a, Vf — 4.436 and «iin ~ 3.330; for 
noise acting on /3 the linear speed is the same while Vf — 4.406; 
for noise acting on 7, uun — 1.235, while vi — 3.075. 



To complete the analysis related to the noisy dynam- 
ics we have also considered Map NB with noise on the 
parameter 7 only. In the pulled case, similarly to the 
Map NA case, we observe saturation of the root mean 
square displacement A(t). In the pushed case A(t) in- 
deed grows in time, but within our time window and for 
the examined parameters {a — 1.9, (3 = 2.5, 70 = —0.5 
and A = 0.4, with Wf = v^n = 3.2046), we do not observe 
any clear scaling. 

Finally we address the chaotic case. Our results are 
reported in Fig. O. 



ficiently large diffusive couplings. Therefore in our case 
we always expect that for chaotic maps A(i) cannot grow 
indefinitely in time, and saturation is the only behavior 
observable. This indicates that the chaotic evolution can- 
not be simply reduced to some erratic behaviour sharing 
common dynamical features with the noisy systems. 



V. CONCLUDING REMARKS 




t 

FIG. 11. Root mean square displacement A(t) as a func- 
tion of time for Map CB with u^, = 1. The values of the 
parameters are in this case a — 1.1 (pushed) and a = 1.8 
(pulled), 7 = 2.5, 13 = -2, and D = 4. The length of the 
chain was L — 60, 000 and the average was performed over 
2,000 initial conditions. 

In this situation, in contrast with the stochastic case, 
we observe for the pulled as well as for the pushed regime 
that A(i) saturates to a constant values. We have verified 
this for Map C5 with W3 = 1 (see Fig. O). An analogous 
result has been observed by Rudzick et al. [|lj for a Cou- 
pled Map Lattice (CML) model. In that case a diffusive 
behaviour for the wandering of the front position is ob- 
served for small diffusive coupling, while for sufficiently 
strong coupling A(i) is shown to saturate. The authors 
of ||l^ argue that this result can be explained as follows. 
For small couplings, the spatial correlation between ad- 
jacent sites is negligible on a distance corresponding to 
their lattice spacing. Therefore, the successive positions 
occupied by the front in its time evolution will be com- 
pletely decorrelated, and the front dynamics induced by 
the local chaotic evolution can be described in term of 
a stochastic process. In contrast, when the coupling be- 
comes large enough, the lattice sites become correlated 
on a length scale larger than the lattice spacing, and this 
implies that a description of the front evolution in terms 
of a stochastic process in not appropriate any longer, not 
even for chaotic maps. Notice that the CML model is a 
spatially discrete model, which can be considered a fair 
approximation of a continuous space system only for suf- 



In this paper we have studied the front propagation 
in reaction diffusion systems with a periodically forced 
reaction term. For this case we have been able to rewrite 
the evolution of the system as a convolution of a spatial 
propagator and a discrete-time map. By implementing 
a quite efficient algorithm for the evaluation of the con- 
volution integral, we have numerically studied the front 
dynamics for deterministic non-chaotic and chaotic maps 
as well as for noisy maps. An analogy with the usual 
FKPP problem can be drawn for non-chaotic maps, al- 
lowing to obtain the expression for the lower and upper 
bound for the speed of the front. Moreover, even when 
the reaction dynamics presents an unstable fixed point 
coexisting with a chaotic (or noisy) behaviour, the anal- 
ogy with the FKPP problem still holds, once the chaotic 
(noisy) phase is identified with the stable fixed point of 
the usual FKPP reaction term. In the pulled regime 
the presence of chaos (noise) plays a poor role and the 
front speed Vf is essentially determined by the dynamical 
evolution around the unstable fixed point. On the con- 
trary in the pushed case Vf depends in a nontrivial way 
on the details of the chaotic (noisy) behaviour. Unfor- 
tunately, the relationship between the chaotic properties 
{e.g. the Lyapunov exponent) and the values of the ob- 
served Uf does not appear simple. In particular, the effect 
of chaos seems to be much stronger than that related to 
the noise: e.g. V{ changes noticeably for chaotic reaction 
dynamics at varying the Lyapunov exponent, while it re- 
mains almost constant as a function of the strength of 
the noise. In contrast, the noise induces a wandering of 
the front around its average position, which is diffusive in 
the pushed case and sub-diffusive in the pulled one. Our 
results confirm the same scaling laws found in continuous 
models with multiplicative noise |l^-|l4| for the time evo- 
lution of the mean square displacement A(t) associated 
to the front position fiuctuations, namely A{t) ^ t^^^ in 
the pushed case and A(i) ~ t-^^^ in the pulled situation. 
This suggests that such results are universal in that they 
do not depend on the details of the system. When a 
chaotic dynamics is considered both in the pushed and 
in the pulled case a saturation of A(t) to a constant value 
is observed, and this is consistent with previous results 
obtained for coupled map lattice with strong diffusive 
couphngs [0. 
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APPENDIX A: DETAILS ON THE 
INTEGRATION ALGORITHM 

Let us now present the numerical details concerning 
the integration of Eq. (@), 



M 



u{x,t + At) 
1 



y/AnDAt J- 



dy e -iD^Mp^ul^x — y,t)] 



(Al) 



which represents the convolution between the field at 
time t + e {i.e. u{x,t + e) = F[u{x,t)]) and the Gaussian 
Kernel 



K{x) 



y/AnDAt 



g 4DAt 



(A2) 



The integration ( [Al| ) can be performed on a computer 
once the field is discretized on a g rid of resolution Ax. 
The integral reported in Eq. (Al) reduces then to the 
sum 

u{iAx,t + At) = 

N/2 

Y^ K{mAx)F[u{{i-m)Ax,t)] , (A3) 

m=-N/2 

where the system length is L = NAx and periodic 
boundary conditions are assumed along the chain. 

In order to im prov e the integration speed we have re- 
stricted the sum (A3) to a small number M of neighbours 
of the site i and to ensure a good precision the true kernel 
has been substituted by a "modified" kernel C{m), cho- 
sen t o m inimize the integration error. With these choices 
Eq. (A3) now reads: 



u{iAx,t + At) 



M 



2_] C{m)u{(i — m)Ax,t + e) 



(A4) 



m=-M 



The problem is now to determine the coefficients {C(?7i)} 
for m = —M, . . . , M . In order to do this, we first suppose 
that an appropriate Fourier basis made up of {2M + 1) 
elements, able to well approximate the field on the chosen 
grid, can be written as 

{exp[j(fca)x]} , 

where the parameter a will be determined later. By 
rewriting Eq. (A4) on this basis one is left with the 
following set of equations: 



^ C(?Ti) exp[— j(fca)?TiAa;] 

= exp[-(fca)2i:»At] 



-M 



(A5) 



Due to the Kernel symmetry we can reduce the system 
of (2M -I- 1) equations ( |A5| ) to a set of M equations. 
This because the elements of the Kernel are symmetric 
around to = (i.e. C{m) — C{-~m))^ and are related by 
the normalization condition 



C(0) 



1 






We can determine the M independ ent elements of the 
Kernel C{m) by solving the system (A5) as a function of 
the parameter a, once the integration time step Ai, the 
diffusion coefficient D, and the spatial resolution Ax are 
fixed. The optimal choice of the parameter a is achieved 
by requiring that the first six cumulants of the discretized 
Kernel reproduce those of the true Kernel within a pre- 
cision of one part over a million and that the quadratic 
sum X]m=-M^^("^) is minimal. In the present paper, 
we have always used D — A, At = Ax = 1 with M — 15 
and a parameter value a — 0.13. 

This integration scheme has been previously intro- 
duced as a possible alternative to the usual pseudo- 
spectral codes to evaluate the dynamical evolution of the 
complex Ginzburg-Landau equation and of the Fitzhugh- 
Nagumo equation in one and two dimensions pl[|. 



APPENDIX B: THE SPEED BOUNDS FOR 
DISCRETE TIME MAPS 



The linear speed (16) for discrete maps can be obtained 
with simple consideration just following the standard rea- 
soning used for the derivation of wijn in the continuous 
time limit: 



d_ 

di 



u^ DV'^u + f[u) 



(Bl) 



As we have already mentioned, the leading edge of the 
propagating front has an exponential shape: 



m(x, t) 



-fj,x+X{fj,)t 



(B2) 



Inserting (B2) in (Bl) and linearizing around u = 
one obtains: 



A(m) - Df,^ + /'(O) 



(B3) 



A stationary phase argument gives a selection criterion 
which allows for the determination of the front speed as 



xi^l) 



t;u„ = min^^ = 2/Dr(0) 



(B4) 
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Let us consider now the discrete-time reaction case, 

-u^DV\+ J2 9{u)5{t-n) , (B5) 

n— — oo 

where we have adopted At = 1. Indicating with F{u) 
the reacting map one gets: 

u{x,t + 0+)^ F{u{x,t)) , 

and integrating the diffusion equation dtu — DV^u be- 
tween i -I- 0+ and t + I one obtains: 



u(x, t + f ) = 



f 



V'^nD 



e 'io F {u{x - w,t))dw , (B6) 



Assuming the shape (B2) and hnearizing around u — 0, 
i.e.: F{u) ~ F'{0)u, a simple Gaussian integration gives: 

The above result implies 

A(/i) =lni^'(0) + Z?/i^ , 



which is nothing but Eq. (B3) now with ln_F'(0) in place 
of /'(O). The same selection criterion gives 



z;iin = 2V^lnF'(0) . (B7) 

In order to estimate Umax hi the discrete case, let us 



shape of u{x,t) given by (B2). It is straightforward to 
show that 



consider again the continuous equation (Bl) with the 

It is 

fin) 



X(n) < Dn^ + r 

0<M<1 u 

From the above inequality it can be shown that 



(B8) 



2\ D max 



fiu) 



0<M<1 U 

In the discrete case we write the following inequality 



(B9) 



u(x,t+l) < 



1 



F{u) 



X / e '^o ii(^x — w,t)dw 



(BIO) 



and assuming that the leading edge has the shape re- 
ported in (B2) one obtains 



A(/i) < Df? 



In 



F{u) 
max 

0<M<1 u 



(Bll) 



This equation is analogous to (B8) therefore in the dis- 
crete case the upper bound for the speed is now 



2W-Dln 



F{uy 
max 

0<M<1 u 
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